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We study the central (inner few hundred parsecs) stellar populations of four starburst 
galaxies (NGC34, NGC1614, NGC3310 and NGC7714) in the near-infrared (NIR), from 0.8 
to 2.4/im, by fitting combinations of stellar population models of various ages and metallici- 
ties. The NIR spectra of these galaxies feature many absorption lines. For the first time, we fit 
simultaneously as much as 15 absorption features in the NIR. The observed spectra are best 
explained by stellar populations containing a sizable amount ( 20 to 56 % by mass) of 1 Gyr 
old stellar population with Thermally Pulsing-Asymptotic Giant Branch stars. We found that 
the metallicity of the stars which dominates the light is solar. Metallicities substantially differ- 
ent from solar give a worse fit. Though the ages and metallicities we estimate using the NIR 
spectroscopy are in agreement with values from the literature based on the UV/optical, we 
find older ages and a larger age spread. This may be due to the fact that the optical is mostly 
sensitive to the last episode of star formation, while the NIR better maintains the record of 
previous stellar generations. Another interesting result is that the reddening estimated from 
the whole NIR spectrum is considerably lower than that based on emission lines. Finally, we 
find a good agreement of the free emission line spectrum with photoionization models, using 
as input spectral energy distribution the synthetic composite template we derived as best-fit. 
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1 INTRODUCTION 

Star-forming galaxies (starburst/H II) are among the best laborato- 
ries to study the evolution of massive stars as well as the physi- 
cal processes that are associated with the earliest stages of galaxy 
fo rmation. The study of this kind of object starts wi th the work 
of lPastorizal ( Il967lll975t) and lSargent & Searld ( Il97(]| ). where it is 
suggested that some galaxies had their (nuclear) spectrum similar 
to those of H II regions. Star-forming galaxies are easily recognized 
from their preeminent emission line spectra, dominated by hydro- 
gen and helium recombination lines as well as strong forbidden 
lines of sulfur, oxygen and nitrogen, among others. These sources 
have also a prominent absorption line spectrum. 

The Equivalent Widths (Ws) of the absorption features give 
insights on the stellar populations of the host galaxy, their study 
is a critical step in the analysis of the continuum and emission 
lines present in the spectra of galaxies. Through the analysis of 
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the stellar content we get information about critical processes 
such as recent episodes of star formation and the evolution- 
ary history of the galaxy. The stellar populations of starburst 
galaxies and the central regions of Seyfert galaxie s have been 
mostly studie d in the UV and optical bands (e.g. Bica l 1 19881: 

Schmitt et al] [l996: Gonzalez-Delgado et al. 1998; Raima nn et al.l 

I II I 1 " .... 
2000alb; Bonatto et all |2000|: ICid Fernandes. Leap & Lacerdd 

| 2003i: iGonzalez-Delgado et al.l |2004| : IWestera et alj |2004 

ICid Fernandes et alj2005 , and references therein). 

The first use of th e near infrared (NIR) dates back to 
nearly thirty years ago. Rieke et al. I ( IT98Q) employ NIR spec- 
troscopy to study NGC253 and M82 and report the detection 
of strong CO band in these sources, which suggests the pres- 
ence of a large population of red giants and supergiants in the 
nuclear region of both sources. Since the, other groups have 
used the ne ar-IR to study star formation, in most cases the CO 
bands (e.g. [Origlia, Moorwood & Clival Il993l : I01ivaetaljri995l: 
lEngelbrac ht et al.ll998l : lLan con et al.'200l'. and references t herein) 
or photometric methods (e.g. Moorwood & Glass 1982; Hu nt et al] 
|2003|) . but overall the whole NIR stellar spectrum (0.8 to 2.4/im) 
with its many absorption features remains poorly explored. 

The obvious reason to use the NIR for starburst galaxies is 
that this spectral range is the most suitable one to unveil the stellar 
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populations in liighly obscur ed sources. However, tracing star for- 
mation in tlie NIR is difficult jOriglia & Olival2000h . To our knowl- 
edge, in this spectral region, except for a few indicators such as the 
met hods based on the CO(2-0) first overtone or the Br7 emission 
(e.g. lOriglia, Moorwood & Olivall 19931 : lOliva et alJI 1995h . the de- 
tection of spectral features allowing the identification and dating 
of stellar populations, mostly those residing in the inner few tens of 
parsecs of galaxies, remains an open question. As stellar absorption 
features in the infrared are widely believed t o provide means for 
recognizing red supergiants jOliva et aEll 19951) . they arise as prime 
indicators for tracing starburst in galaxies. 

Besides the short-lived red supergiants typical of very young 
stellar populations, the NIR hosts also the contribution from the 
thermally- pulsating asymptotic giant branch (TP-AGB) stellar 
phase, which is enhanced in young to interme diate age stellar pop- 
ulations (0.2 4: t ^ 2 Gvr. Warasto"rilll998h . The contribution of 
this stellar phase in stellar population models has been rec ently in- 
cluded in both the energetics and the spectral features by (Maraston 
[2005a hereafter M05). In particular, these models employ empirica l 
spectra of oxygen-rich and carbon stars ( lLancon&Woo3l2000h . 
which are able to foresee characteristic NIR absorption features 
such as the 1.1 /im CN band. TP-AGB stars leave a unique fin- 
gerprint on the integrated spectra (cfr. M05, Fig. 14), hence when 
detected they can help the age dating of stell ar populations from 
integrated light (M05). In fact, very recently iRiffel et al] ^20071) 
detected the I.I/^m CN band in the spectra of active galactic nu- 
clei and starburst galaxies. In addition, |Dottori e t al. (2005) propose 
that carbon stars are natural candidates to explain the colour excess 
of the nucleus with respect to the circumnuclear star formation ring 
of NGC 1241. Furthermore, the typical NIR spectral shape due to 
the presence of TP-AGB stars has been detected in Spitzer data for 
high-redshifted galaxies ( Maraston et al. 2006) . 

The knowledge of the stellar populations that dominate the 
light in the NIR becomes very important nowadays as the spec- 
troscopic observation of galaxies is being directed to this spec- 
tral region. With the new generation of infrared arrays, it is now 
possible to obtain spectra at moderate resolution on faint and ex- 
tended sources. Also, the availability of cross-dispersed spectro- 
graphs offering simultaneous wavelength coverage in the interval 
0.8—2.4 /im, allow the study of the NIR region avoiding the aper- 
ture and seeing effects that usually affect JHK spectroscopy done 
in long-slit mode and single-band observations. 

With this in mind, we present for the first time a detailed fitting 
of the whole 0.8-2.4/im spectral range, for the Starburst/H II galax- 
ies studied by Riffel, Rodriguez- Ardila & Pastoriza ( 2006, here- 
after RRP06). This paper is structured as follows: The data are pre- 
sented in Sec. [2] In Sec. [3] we describe the fitting method and how 
the measurements are done. Results are presented and discussed in 
Sec.|4] In Sec.|5]we focus on the free emission line spectrum. The 
final remarks are given in Sec.|6] 



2 THE DATA 

For this work we choose NGC 34, NGC 1614, NGC 3310 and 
NGC 7714, which are typical Starburst/H II galaxies and are 
very well studied in other spectral ranges. The NIR spectra of 
these sources, presented in RRP06, were obtained at the NASA 

3 m Infrared Telesc ope Facility (IRTF). The SpeX spectrograph 
dRavner et al]|2003h . was used in the short cross-dispersed mode 
(SXD, 0.8-2.4/im). The detector used was a 1024x1024 AL- 
ADDIN 3 InSb array with a spacial scale of O.I5'7pixel. A 



0.8"xI5'" slit was employed giving a spectral resolution of 360 
kms^^. For more details about instrumental configuration and 
sample selection see RRP06. A rapid inspection of Fig. 13 in 
RRP06 allows us to the see that all the spectra are dominated by 
strong absorption features due to stars. For reference, we plot in 
Fig.[T]the most proeminent ones observed in NGC 7714, including 
the molecular bands of CN, CO and those of atomic transitions of 
Ca II, Ca I, Na I, Si I and Mg I. Emission lines present in the NIR 
spectrum of this object are also identified. 



3 ANALYSIS OF THE STELLAR POPULATIONS : 
METHOD AND MEASURES 

In order to describe the age distributions of the stellar populations 
of the sources and their metallicities, we compare the observed 
spectra with various combinations of stellar population templates 
of various ages and metallicities. In th e ana lysis we use a method 
similar to the one introduced bv iBical jl988 ). The latter employs 
empirical star cluster spectra and combine them in order to estimate 
the combination of ages and metallicities that allows the best match 
to the galaxy spectra. Here instead, we use as input for the base 
of elements the synthetic Simple Stellar Population (SSP) models 
computed by M05. 

The fitting technique used here exploits the equivalent widths 
(W\) of spectral absorption features and the continuum fluxes {Fx) 
at selected wavelengths intervals, free from emission/absorption 
lines. These quantities are measured in the galaxy spectra and com- 
pared with combinations of the base of elements, which in turn have 
known ages and metallicities. An algorithm generate these combi- 
nations and compare the composite synthetic W\ and -Fa with the 
observed values. In practice, the algorithm tries to find the mini- 
mum of the function, 

/(X)= ^[Wgal(i) - Wsp(i)]' + £[Fgal(j) - Fsp(j)]' (1) 
i=l j=l 

were Wgai and Wgp are, respectively, the equivalent width of a 
given line / of the galaxy and of a Stellar Population model (here- 
after SP). The Fgai and Fsp are the continuum fluxes at a selected 
point i of the galaxy and of a SP, respectively. The possible solu- 
tions are those that reproduce, within allowed limits, the W\ and 
Fx observed in the galaxy spectra. The results of this matching pro- 
cedure are flux fractions at a selected A of the different components 
of the base of elements. 

In the present work we have upgraded the algorithm of 
ISchmitt et al.l ( 199^ in order to include a search for the internal 
reddening of the stellar population E(B-V)sp for F\ points, free 
from emission/absorption lines, located in the NIR region. The 
ICardelli. Clayton & Mathid ( Il989[) [CCM] reddening law was em- 
ployed for this purpose. 

We have defined the NIR continuum fluxes at 0.81/im, 
0.88/^m, 0.99/im, I.06^im, 1.22^im, 1.52/im, I.70^im, 2.09/im, and 
2.19/im, all normalized to 1.22/im, to be used as input parameters 
for the algorithm. As can be observed in Fig [2] we have used an 
average value of the continuum points located between the verti- 
cal lines. The measured values in the galaxy spectra as well as the 
contin uum intervals are presented in Tab[T] The PACCE iVale et all 
I2OO8I) code was used to measure the fluxes. 

Bandpasses and continuum points used to compute the W\ 
of the strongest absorption lines were taken from the literature or 
newly defined here. For the CaT lines and for the CN band we 
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Figure 1. The most prominent absorption and emission lines observed in the spectrum of NGC7714. Note that here the spectral resolution was degraded in 
order to match the M05 models (see Section[3j. 



Table 1. Defined continuum fluxes, normalized to unity at 1.2230/xm, mea- 
sured in our galaxy sample. The errors on _F;^ are 2% in all cases. 



Interval (^tm) NGC34 NGC 1614 NGC3310 NGC7714 



0.813-0.816 


0.86 


0.76 


1.29 


1.41 


0.879-0.881 


0.94 


0.88 


1.27 


1.41 


0.993-0.995 


1.02 


1.00 


1.19 


1.25 


1.057-1.060 


1.03 


1.03 


1.13 


1.18 


1.220-1.226 


1.00 


1.00 


1.00 


1.00 


1.515-1.525 


0.91 


0.85 


0.79 


0.77 


1.699-1.703 


0.92 


0.84 


0.72 


0.69 


2.085-2.100 


0.63 


0.54 


0.41 


0.39 


2.185-2.195 


0.56 


0.49 


0.37 


0.34 



and 1.23 fim. The continuum was also adjusted as a cubic spline 
for the measurements of Mg I, Mg 1 1 .577 1 ^m, Si 1 1 .5894 ^m 
and CO 1.6175 /im using points between 1.43 fim and 1.69 fim. In 
the same way, we have adjusted the continuum as a cubic spline 
for the lines located in the if -band, using continuum points free 
from contaminations of emission/absorption features in the 2.09- 
2.37 /im spectral range. In Fig. [3] we illustrate the continuum and 
bandpasses used to compute the W\. 

Note that the resolution of the observed spectra has been de- 
graded (see Fig.[T]( to that of the M05 model This procedure was 
done using the task gauss of the IRAlfl software. The measured 
values, as well as their uncertainties are listed in Tab.|2] Errors were 
estimated using the standard deviation of three different choices of 



emplo yed the definitions of iBica & AUoinl h987h and lRiffel et al.l 
( l2007l) . respectively. The W\ were measured using the band- 
pass as Hsted in Tab. |2] For Alll.1285 /im, Nall.l383/im and 
Si 1 1.21 12 /im the continuum was adjusted as a cubic spline us- 
ing points free from emission/absorption lines between 1.02 /im 



^ 20 A at the bluest lambdas of our spectram up to 1 fim, 50 A from 1 fim 

to 1.6 fim and 100 A up to the end of the galaxy spectrum. 

^ IRAF is distributed by the National Optical Astronomy Observatories, 

which are operated by the Association of Universities for Research in 

Astronomy, Inc., under cooperative agreement with the National Science 

Foundation. 
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Figure 2. Continuum intervals free from emission/absorption lines of 
NGC 7714. We use average values of the points between the vertical lines. 




2.2 2.26 



A(/im) 



Figure 3. Bandpasses (dashed lines) used to compute the equivalent widths 
of NGC 77 1 4. Al I 1 . 1 250^m and Si I 1 .2 1 1 2^lm are from NGC 34. The red 
line (full line) is the ajusted continuum. 



the continuum level around a given W\, which basically reflects 
the subjectiveness associated with continuum determinations. 

The spectral base was taken from the evolutionary population 
synthesis (EPS) models of M05. These models are particularly in- 
teresting for studying the stellar populations in the NIR because, 
as discussed in the Introduction, they include empirical spectra of 
carbon and oxygen rich stars evolving through the TP-AGB phase. 
Thus, the models can predict the stre ngth of molecular features like 
CH, CN and C2. iRiffel et all ( |2007|) report, for the first time, the 
detection of the 1 . 1 /xm CN band in a sample of galaxies including 
the sources studied in the present work. CN is particularly strong 



in young/intermediate stellar populations with ages in the range 
^0.3^ t ^2 Gyr, though is also detected in red supergiant s of mas- 
sive stars (e.g. lLancon & Wooj|2000| ; lLancon et al.ll2007h . 

As base set we consider solar metallicity SSPs with a Salpeter 
initial mass function and the following age distributions: 0.001, 
0.01, 0.03, 0.05 0.2, 0.5, 0.7, 1 and 13 Gyr The ages are chosen 
taking into account the following: (i) very young starbursts (1 Myr 
population) can be present; (i/) re d supergiants (RSG) are impor- 
tant c ontributors to the NIR (e.g. lOliva et al] 1 19951 : iLancon et al.l 
lloml and are included in the 10, 30 and 50 Myr populations fol- 
lowing the input stellar tracks of the M05 models (from the Geneva 
databa se, see M05 for re ferences) ); (Hi) stars in the TP-AGB phase 
(M05, lRiffel etal.ll200l populations with 200, 500, 700 Myr and 
1 Gyr). We additionally include an old model 13 Gyr The behavior 
of the Wx with age is presented in Fig. |4] The base set is capa- 
ble to discriminate between the stellar populations in the RSG and 
TP-AGB phases (see Fig.|4]l. Here it should be noted that while the 
M05 models include empirical spectra of C-O-rich TP-AGB stars, 
they do not include empirical spectra of re d supergiants. The la tter 
show a large scatter in observed properties jLancon et alj2007h . In 
addition, also the energetics of the RSG phase is not well known as 
it depends on mass-loss in massive stars, which is in turn highly un- 
certain. On the other hand, RSG display near-IR absorption bands 
and a proper spectrum should be used. To explore how much we 
could be in error by neglecting the true spectrum of a RSG, we 
made a 30 Myr, M05 model in which we include a stellar spectrum 
of an 0-rich star from the database of lLancon & Wood ( 2000), that 
displays very similar features to a RSG (see lLancon & Woodl200d 
Fig. 7). We have then repeated the fitting procedure for the galaxy 
NGC 77 14 using this 30 Myr model instead of the standard one. 
We find that the results are essentially the same. The measured W\ 
and F\ for the SSPs used in the base are given in Tab.[3]and Tab|4] 

It is important to note that we have tested other metallicities 
in the range Zq Z ^ 2^©. To find the best solution we 
compute the x of the residuals of the W\ and F\ values, i.e. we 
solve the equations: 



(2) 



and 



Nf 



(3) 



The and AFa are the differences between the observed 

and synthetic values for the W\ and Fx-, respectively. As can be 
seen in Tab. |5] we find that the solar metallicity provides the best 
solution for the four galaxies. This result is in agree ment with the 
analysis of the optical re gion of spiral bulges (e.g. ISchmitt et al.l 
ll99d : ISaraivaetalJl200lh . Hence, as a first result of this paper, we 
find that the nuclear regions of these starburst galaxies are charac- 
terized by a solar-like metallicity. 



4 RESULTS 

The results of our fitting procedure, in terms of flux fractions at 
1.2230 /im, are shown in Fig [5] and Tab. [6] As can be seen, there 
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Table 2. Equivalent widths (in A) of the absorption bands measured on the sample galaxies. 



Feature^ 


A(^m) 


Band Pass (pm) 


NGC 34 


NGC 1614 


NGC 3310 


NGC 7714 


CaTi 


0.8498 


0.8476-0.8520'" 


3.75±0.06 


4.57±0.07 


2.44±0.21 


2.36±0.02 


CaTz 


0.8542 


0.8520-0.8564*" 


3.89±0.01 


4.75±0.05 


3.87±0.18 


3.40±0.03 


CaTg 


0.8670 


0.8640-0.8700*" 


3.70±0.16 


4.05±0.09 


3.12±0.16 


3.31±0.13 


CN<= 


1.0950 


1.0780-1.1120 


24.72±0.13 


25.77±0.17 


16.00±0.37 


24.34±0.16 


All 


1.1250 


1.1200-1.1300 


0.45±0.04 








Nal 


1.1395 


1.1335-1.1455 


1.55±0.10 


2.76±0.08 


3.10±0.11 


1.16±0.17 


Si I 


1.2112 


1.2025-1.2200 


3.57±0.50 








Mgl 


1.5771 


1.5720-1.5830 


5.64±0.30 


6.34±0.40 


5.05±0.09 


4.42±0.04 


Si I 


1.5894 


1.5870-1.5940 


1.93±0.03 




3.82±0.08 


2.01±0.07 


CO 


1.6175 


1.6110-1.6285 


4.27±0.40 


5.37±0.15 


4. 85 ±0.24 


4.15±0.08 


Nal 


2.2063 


2.1936-2.2150 


4.78±0.05 


3.95±0.14 


2.87±0.12 


2.44±0.06 


Cal 


2.2655 


2.2570-2.2740 


4.48±0.04 


4.07±0.02 


2.36±0.13 


3.57±0.01 


CO 


2.2980 


2.2860-2.3100 


16.72±0.29 


17.66±0.45 


13.87±0.44 


14.22±0.32 


CO 


2.3255 


2.3150-2.3360 


8.91±0.14 


9.75±0.16 


9.06±0.35 


9.24±0.13 


CO 


2.3545 


2.3420-2.3670 


14.68±0.50 


16.02±0.15 


15.74±0.12 


13.20±0.04 



a) The identifications of the atomic hnes are based on 
Gushing. Ravner & Vacca (2005) . 

b) From Bica & Alloin (1987). 

c) For details see Ritfel et al. (2007) . 



Table 3. Equivalent widths (in A) for the SSPs used in the base. 



Table 5. x values for the different metaUicities. 



Age"- CuTi CaT2 CaTs CN 



CO Nai Cai CO CO 



0.001 1.45 1.49 1.88 0.00 0.00 0.00 0.00 0.00 0.30 0.00 1.29 0.00 0.00 0.00 0.00 

0.01 8.70 7.82 4.53 0.00 0.00 1.15 0.00 9.24 5.12 0.74 2.81 0.00 20.66 13.81 18.02 

0.03 2.70 2.81 3.35 12.52 0.00 0.4! 0.00 0.25 1.81 1.06 1.60 0.00 7.53 6.13 6.33 

0.05 2.56 2.75 3.20 11.26 0.00 0.57 0.00 0.62 1.81 0.84 1.59 0.00 7.79 6.20 6.07 

0.2 2.33 2.59 3.00 10.84 0.00 0.71 0.00 1.10 2.0! 0.68 1.62 0.00 8.30 6.59 5.25 

0.5 3.08 3.38 2.16 31.50 0.00 2.25 3.59 4.91 0.77 2.98 2.88 4.91 18.5! 11.97 12.69 

0.7 2.67 3.00 1.64 34.63 0.00 2.89 4.3! 5.77 1.49 3.76 3.47 6.70 21.10 12.33 12.22 

1 3.16 3.48 1.91 34.65 1.64 3.07 4.74 6.86 2.26 4.40 4.58 7.55 22.13 10,54 10.87 

13 4.07 3.76 2.78 0.00 0.00 1.24 0.22 3.80 1.92 0.63 2.12 0.00 13.71 9.51 7.10 



Source 



2Z. 



© ^0 



^0 



50 



2Z. 



X(Fx) 



^0 



^0 



NGC 34 2.26 

NGC 1614 2.90 

NGC 3310 2.22 

NGC 7714 2.13 



1.08 
1.53 
1.54 
1.35 



1.39 
2.28 
1.72 
1.40 



1.69 
2.14 
2.31 
5.82 



0.21 
0.11 
0.06 
0.03 



50- 



^0 



0.10 0.22 0.23 

0.03 0.11 0.18 

0.04 0.08 0.07 

0.02 0.04 0.04 



Table 4. Continuum fluxes, normalized to unity at 1.2230^m, measured on 
the base SSPs. 



Age 


0.81 


0.88 


0.99 


1.06 


1.22 


1.52 


1.70 


2.09 


2.19 


0.001 


4.57 


3.48 


2.21 


1.74 


1.00 


0.43 


0.28 


0.28 


0.10 


0.01 


1.40 


1.50 


1.49 


1.32 


1.00 


0.81 


0.76 


0.76 


0.35 


0.03 


2.57 


2.36 


1.75 


1.47 


1.00 


0.57 


0.44 


0.21 


0.18 


0.05 


2.42 


2.18 


1.67 


143 


1.00 


0.60 


0.46 


0.23 


0.19 


0.2 


2.21 


1.98 


1.58 


1.38 


1.00 


0.62 


0.47 


0.47 


0.20 


0.5 


1.50 


1.51 


1.28 


1.24 


1.00 


0.69 


0.68 


0.68 


0.33 


0.7 


1.47 


1.49 


1.28 


1.24 


1.00 


0.70 


0.71 


0.71 


0.35 


1 


1.45 


1.49 


1.29 


1.24 


1.00 


0.71 


0.72 


0.72 


0.36 


13 


1.58 


1.50 


1.37 


1.26 


1.00 


0.74 


0.62 


0.62 


0.28 



is a good agreement between the observed spectrum and the syn- 
thetic one for all galaxies. We also list the derived E(B-V) as well 
as the morphological type of the sources, as obtained from NED. 
The stellar population in terms of mass fraction was calculated us- 
ing the mass-to-light ratio (M/L) in the J-band[f|, assuming that 
the J-band (centred at ~1.2/im) M/L is equal to M/Li. 2230^111, 
where Li, 2230pm is the continuum luminosity at 1.2230/im. The 



see www.dsg.port.ac.uk/ marastonc/ 



mass fractions are listed in Tab. [6] We compare the flux and mass 
fractions of each population for the four sources in Fig. |6] In the 
remaining of this section we comment on the results for each indi- 
vidual source. 

4.1 NGC 34 

The optical properties of NGC 34 (Mrk93 8, MCG-02-0 1-032) 
are di scussed in detail in the recent work of ISchweizer & Seitzej 
toot . Using photometry and Lick indices (when spec- 
troscopy was poss i ble) a nd comparing them with the models of 
iBruzual & Charloll ( l2003h they conclude, that NGC 34 has a young 
stellar disk and supports a rich system of young star clusters. In 
fact, as can be observed in Figs. |5] and [6] and in Tab. [6] the light at 
1 .2230 /im in the 230 central pc of NGC 34 is also dominated by a 
stellar population with solar metallicity and a young-intermediate 
aged population (~lGyr). These findings are in full agreement 
with the fact that NGC 34 has in i ts spectrum a strong 1 . 1 /im CN 
absorption band jRiffel et alj2007h characteristic of stellar popula- 
tions with that age (M05). 

The colour excess E(B-V) of this sou rce was calculated by 
iRodrfguez-Ardila. Riffel & Pastorizal ilOOj. hereafter RRP05) us- 
ing the emission line ratio Pa/3/Br7. They obtain E(B-V)=1.32 
mag. As can be observed, the E(B-V) estimated by means of the 
emission lines is larger than that obtained from the spectral fitting 
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Figure 4. Absorption Unes as functions of age. Note tiiat tiie 2.3/xm CO bands as well as the CN are plotted with a different scale for Wx. 



Table 6. Results in terms of flux and mass fractions (first and second Une, respectively). 



Source Morphology^ E(B-V) 



% of age contributions ^ . 



(mag) 0.001 0.01 0.03 0.05 0.2 0.5 0.7 1 



13 



NGC34 Sc 0.95±0.05 1±2 5±5 8±7 7±7 7±7 6±7 10±9 53±7 3±2 

1±4 Oil Oil Oil 7±14 4±10 8±15 56±14 23±31 

NGC1614 SB(s)c 1.00±0.06 1±1 1±4 6±2 12±3 6±3 4±8 6±5 43±7 21±8 

Oil ()±0 0±0 0±0 3±3 1±7 2±5 20±9 73±72 

NGC3310 SAB(r)bc 0.55±0.03 2±1 12±6 13±10 Bill 12±11 4±5 6±7 30±6 8±4 

1±1 1±0 1±0 1±0 10±16 2±5 4±9 27±9 54±45 

NGC7714 SB(s)b 0.46it0.02 7±3 7±4 10±8 13itl0 11±9 4±4 6it5 34±4 8±3 

3±6 0±0 0±0 0±0 9±16 2±5 4±8 30±8 51±42 



t From NASA/IPAC Exttagalactic Database. 
% Ages in Gyr. 
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(Tab.[6l(. The difference can be due to the fact that Pa/3 is strongly 
affected by absorption lines like Til 1.2827/im and Til 1.2851/im 
as well as Pa/3 itself, as observed in Fig.|5] The discrepancy can 
be also explained considering that hot ionizing stars can be associ- 
ated to a dustier region with respect to the c old stellar population 
jCalzetti. Kinney & Storchi-Bergmannlll994l) . 



4.2 NGC 1614 

This source (Arp 186, Mrk 617,11 Zw 015) is considered as a labora- 
tory f or studying the evolution of a starburst jAlonso-Herrero et alj 
l200lh . The HSTfNlR camera a nd multiobject spectrometer (NIC- 
MOS) observations reported bv lAlonso- Herreroetal] ( l200ll) show 
deep CO stellar absorption, tracing a starburst nucleus with a di- 
ameter of ~45 pc, surrounded by a ~600 pc diameter ring of ex- 
tremely luminous supergiant Hll regions. The stellar population 
that dominates the light in the inner ~9kpc of this sou rce at 
4020A is studied by ICid Fernandes. Leao & Lacerdal ( l2003l here- 
after FLL03), using abase of elements compos ed by 12 star cluster s 
of different ages and metallicities, defined by Schmidt et al.l ( ll99lh . 
FLL03 find that the dominant stellar population is 100 Myr old (~ 
39 %). Here we find that the light at 1.2230^m in the inner 154 pc 
is dominated by a stellar population of 1 Gyr age (see Tab.|6]and in 
Fig.|6). The dominant 1 Gyr stellar population found in this source 
as well as in NGC 34 agree with the fact that the two spectra are 
very similar. 

The difference between our results and those of FLL03 is 
probably related to the fact that the aperture used by FLL03 in- 
clude the light of the younger bulge/arm stars. Interestingly, we 
tried a combination of SSPs models in the proportions determined 
by FLL03 in the optical region, but the resulting synthetic template 
did not reproduce the observed NIR spectra of the inner 154 pc of 
NG C 1614. 

IPuxlev & Brand! ( Il994h study the dust extinction of this source 
by means of hydrogen line ratios. They argue that the NIR line ra- 
tios are inconsistent with a foreground screen of dust obscuring a 
compact radiation source, implying a non-uniform distribution of 
dust. They also show that if the dust distribution is homogeneous, 
then the total extinction is well fitted with an A v=15±2.5 mag. 
Using only Pa/3/Br7, with a point source model, IPuxlev & Brand 
( 1 19941) determine E(B-V)=1.56 (using Rt/=3.0). This E(B-V) value 
is ~ 2 times larger than that estimated by RRP05. This galaxy is 
the only one of our sample where the reddening determined by the 
stellar population fitting is smaller than that determined by means 
of emission lines. This testifies the complex nature of the dust dis- 
tribution in this object. This source should be observed with 2D 
spectroscopy in order to map the dust distribution along the galaxy. 

4.3 NGC 3310 

NGC 3310 (Arp 217) is a nearby disk galaxy undergoing intense 
starburst activity. It has been the subject of many studi es across 
the electromagnetic spectrum (see lElinegreen et al.ii2002l . and ref- 
erences therein). This galaxy has a circumnuclear ring of star for- 
mation with angular diame ter of 713 to 770 pc (a ssuming a diame- 
ter between 8" to 12" from lElmegreen et al.l2002l and a distance of 
13.24 M pc from NED). NGC 33 10 hosts a spectacular starburst, as 
seen in numerous H II regions. However, it does not show evidence 
for even weak 0-type st ellar populations as s een from the HUT 
spectrum of this source jLeitherer et al.ll2002h . The nucleus and 
six surrounding H II regions, four of which are located at less than 



400 pc from the nucle us, were studied i n the s pectral range between 
3600A and 9600A bv lPastoriza et al.l ( Il993h . Here we present the 
contributions of the different stellar populations at 1 .2230/im in the 
inner 56 pc of this object. As can be observed in Fig.|6]and in Tab.|6l 
NGC 3310 displays an intense starburst activity, with four bursts 
of star formation. One dominant at 1 Gyr, which contributes with 
30±6% and three remarkable bursts at 30Myr (13±11%), 50 Myr 
(13±11%) and200Myr(12±ll). 

Regarding the reddening, 

iRodrfguez-Ardila. Riffel & Pastorizal ( l2005h has determined 
an E(B-V)=0.76 using emission line ratios. As in NGC 34, the 
color excess determined from the stellar population fitting is 
smaller than that derived from the emission lines of the gas. This 
discrepancy is probably due to the fact that the Pa/3 is affected by 
absorption lines and/or by the dust distribution, as discussed in the 
case of NGC 34. 

4.4 NGC 7714 

NGC 77 14 (Mrk 538 , Arp 284) has been classified by 
IWeedman et al.l imih as a prototypical starburst galaxy. The 
optical stellar populations in the inner ~5 kpc of this source have 
also been studied by FLL03. They find a 100 Myr old (~ 33 %) 
dominant population. As for NGC 1614, we have combined the 
SSPs in the proportion calculated by FLL03 but the synthetic 
template did not match the observed NIR spectrum in the inner 
115pc. 

Our synthesis results indicates that NGC 77 14 is an active 
starburst galaxy. As for NGC 3310, three bursts of star formation 
were found, a dominant 1 Gyr (34±4%) old one and two minor 
bursts with ages 30 Myr (10±8%) and 50 Myr (13±10%). The 
results agree with the facts that this galaxy displays a proeminent 
lAjim CN band (Riffel et al. 2007) and that this source is in merg- 
ing process due to a recent off-cen ter collision (100-200 Myr) with 
NGC 77 15 dStruck & Sm"ifl]|2003l) . 

With very little silicate absorption and a temperature of 
the hottest dust component of 340 K, NGC 7714 is defined by 
ffirandl et al. (2004) as the perfect template for a young, unobscured 
starburst. However, we measured an E(B-V)=0.47 from our spec- 
trum, based on hydrogen line ratios (RRP05). These values agree 
very well with that determined by IPuxlev & Brand! ( Il994t) . The 
E(B-V) determined using the stellar population fitting is similar to 
that determined with the emission lines. 

In concluding this Section, we call the attention to the fact 
that our results are not driven by the CN band. In order to address 
this point we perform the synthesis for NGC 7714 excluding the 
CN band from our fitting. The age distribution found in this case is 
very similar to the one obtained using all absorption lines 0. 



5 IONIZED GAS 

After subtracting the stellar population contribution we are left with 
a pure emission line spectrum allowing the study of the ionized 
gas. We consider that the emission line spectrum of the galax- 
ies is produced solely by photoionization due to stars. This as- 
sumption is based on X-ray observations reported for the four 
sources, which were found to be too weak to be originated by 

4 Age (Gyr) 0.001 0.01 0.03 0.05 0.2 0.5 0.7 1 13 
Cont. (%) 7±3 6±4 4±4 14±9 13±10 3±4 6±6 34±5 11±8 
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Figure 5. Fitting results in terms of flux contribution at 1 .2230/im. The galaxies are identified. Top panel: solid lines represent the observed spectra, shifted 
by a constant and smoothed to the models resolution. Dashed lines represent the sum of the different stellar populations contributions. Dotted Unes show the 
component stellar populations, according to Tab.|6] Bottom panel: the synthetic template is subtracted from the galaxy observed spectrum. 



a buried AGN. In NGC7714, for instance, the X-ray emission is 
expla ined by the presence of supernova remnants ( Weedman et'al] 
Il98lh . The ROSATASCA observations of NGC 33 10 suggests that 
the soft emission is probably a super wind while the nature of the 
hard emission is more uncertain with the likely origins being X- 
ray binaries, inverse Compton scattering of infrared photons, an 
AGN or a very hot ga s component (Zezasetal. 1998). However, 
using ROSAT High Resolution Imager re- 
port that the ultraluminous X-ray sources of this galaxy are located 
around knots on spiral arms. X-ray spectroscopy of N GC 1614 sug- 
gests that this galaxy may harbor an obscured AGN iRisaliti et all 



( I2OOOI) , but VLBI studies with a sensitivity limit of 0.9 mjy do not 
detect a compact radio core in NGC 1614 Hill et al. (2001). In ad- 
dition, the three galaxies listed above are located, in a diagnostic 
diagram envolving NI R emission line ratios, in the region corre- 
spondig to SB galaxies jRodriguez-Ardila et al .|2004 RRP05). For 
the last galaxy, NGC 34, its nuclear spectrum is presently thought 
to be dominated by a highly obscured starburs t, with a likely weak 
AGN contribution JSchweizer & Seitzeill2007h . 

With the above in mind we compute a grid of models of Hll 
regions in order to reprod uce the observed emissi on line intensities. 
Following the method o ors & Copettil ( |2003) . we use the pho- 
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Figure 6. Percentage contributions of each SSP to the total stellar population of the galaxies. Solid lines represent the mass contribution whereas dashed lines 
the flux fraction. Note that the region between 12800 and 13200 Myr is shown in a separate box (right panel) for displaying purposes, as well as that this figure 
is equivalent to Tab|6] but with the ages in Myr 



toionization code CLOUDY/C07.02 (see lFerlanj|2003h . The free 
input nebular parameters of the models are the abundances and ion- 
ization parameter U. The nebula was considered as an expanding 
sphere, with filling factor e=0.01 and a constant electron density 
of Ne = 500cm^^. The spectral energy distribution (SED) used to 
compute the models are the synthetic templates derived from the 
stellar population fitting. We extrapolate the NIR results from the 
ultraviolet to the far-infrared. In this process, we weight the con- 
tribution of each SSPs with the flux contributions listed in Tab. [6] 
and sum them. The final SED is very similar for the four galaxieo 
For this reason, we assume as representative, the resulting template 
obtained for NGC 7714. 

The predicted emission line ratios are in good agreement with 
the observed ones, as can be seen in Tab. |7] Note that we have 
used ([Slll]9531+Pa8)/9069A instead of the [Slll] line ratio be- 
cause the spectral resolution, 20A, in this spectral region does not 
allow deblending [Slll] 9531 A and Pag (A9548 A). In addition, we 
plot in Fig. |7] ([SlIl]9531+Pa8)/Pa/3 versus [Feuj/Bij. The re- 



° The young stellar populations, which ionize the gas, are nearly the same 
for all galaxies. 



suits show that the gas has an ionization parameter in the range 
-3.0 ^ log( U) ^ -2.5, and abundances of S and Fe (X/H) 
between 0.08 and 0.14 solar. The low sulfur a bundances are in 
agreement with the findings o f lKehrig et al.H2006i) . who studied the 
chemical abundances in a sample of 22 Hll galaxies by means of 
the[Siii] 903 1,953 1 A emission lines. They inferred that most of the 
objects in their sample have S/H abundances in the range between 
1/20 and solar. Regarding the [Fell] the ratio [Fell] 1.257 /^m/Br7 
found for our sample varies from 1.40±0. 3I to I.77±0.32 (see 
Tab.|7]and Fig.|7|, similar to t hose reported b y lsimpson et al.l j 1 9961 . 
0.65 and 2.70)11 Models of lColinal ^1993') predict ratios between 
0.13 and 1.90 for starburst galaxies, well within the range of our 
observations. 

It is interesting to note that while the best fitting solutions 
for the stellar populations require SSPs of solar metallicity, the 
S and Fe abundances of the emission gas are lower than solar. 

^ ISimpson et all il99d) indeed reported values of the ratio 
[Fe II] 1.257 /im/Paa in the interval 0.11-0.46 for starburst galax- 
ies. Assuming case B, Paa/Br7=5.88, we estimated that the ratio 
[Fell] 1.257 fim/Paa translates to [Fell] 1.257 /im/Br7 between 0.65 and 
2.70. 



10 Rijfel et al. 



Table 7. Properties of the emission line spectrum. Calculated versus ob- 
served values. 



Feature 


Cal." 


NGC34 


NGC 1614 


NGC 33 10 


NGC 77 1 4 


Pa7/Pa/3 


0.56 




0.53±0.08 


0.61 ±0.05 


0.52±0.03 


Br7/Pa/3 


0.17 




0.16±0.03 


0.19±0.02 


0.18±0.01 


([Sill] 9531+Pag)/9069A 


2.70 


2.71±0.68 


2.70±0.13 


2.67±0.12 


2.48±0.05 


([Sill] 9531+Pag)/Pa/3 


b 


1.41±0.1l 


2.66±0.14 


3.05±0.14 


3.57±0.15 


[Fell]/Br7 


b 




1.40±0.31 


1.77±0.32 


1.47±0.27 



a) Calculated using CLOUDY (see for example ^^^^003). 

b) see Fig.lTl 



One mechanism that can be invoked to explain this difference is 
interaction among galaxies. Three of the four sources observed 
(NGC 34, NGC 1614 and NGC 7714 ) are stro ngly interacting ob - 
jects jThean et alj l2003l : iHattori et al. 2004; Dopita etalj|2002h . 
The gas motions created by the intera ction can significantly al- 
ter the chemical state of the galaxies tKoeppen & Arimotdll99(]| : 
lDalcantorill2007h . and modify the usually smoot h radial metallicity 
gradient often fo und in is olated disk galaxies jHenrv & WorthevI 
[l999). Recently. Ixewlev et al. (2006) found that the O/H abun- 
dance in the central region of nearby galaxy pairs is systemati- 
cally lower than that of isolated objects. These authors suggest 
that the lower metallicity is a consequence of gas infall caused 
by the interaction. In addition, streaming motions of neutral gas 
(Hi) towards th e galaxy center along t he spiral arms was ob- 
served in M 8 1 jAdler & Westpfahj [l996 ) and more recently in- 
fall of mol ecular gas (H2) tow ards the center of NGC 4051 was 
reported by iRiffel et all feOOSh . Similar motions of molecular gas 
towards the nucleus along the bar have been observe d in NGC 4151 
jMundel l & Shone 1999). In the same way. Mazzuca et al.l J2006l) 
have shown that the star-forming ring in NGC 7742 rotates in the 
opposite direction of the bulge stars, and therefore it is made of ex- 
ternal material. Therefore, interactions and inflows of matter from 
the outer low metal abundance medium can dilute the nuclear high 
metallicity interstellar medium. 



6 FINAL REMARKS 

We have analyzed in NIR spectroscopy, from 0.8 to 2.4 /im, the 
stellar populations in the inner few hundred parsecs of the starburst 
galaxies NGC 34, NGC 1614, NGC 3310 and NGC 7714. We use a 
composition of simple stellar population models with various ages, 
from 1 Myr up to 13 Gyr, and we also test the effect of metallicity. 
For the first time, we fit simultaneously as much as 15 absorption 
features in the NIR. We find that all galaxies host a remarkable frac- 
tion of 1 Gyr stellar populations. This result gives support to the 
fact that the 1.1/imCN band is a powerful tracer of inte rmediate 
age stellar populations jRiffel et alj|2007l : lMarastorj|200^ . though 
we have tested that the ages we derive depend on the multiband fit- 
ting rather than on one specific absorption. Our results are in broad 
agreement with similar analysis in the optical region, though we 
tend to find older ages and a wider age spread. At least part of this 
result is due to the fact that the NIR is a region more suitable to pick 
up underlying older stellar populations, and to display the unique 
absorptions proper to 1 Gyr-old stellar populations featured by TP- 
AGB stars. Our finding that the central regions of the SB galax- 
ies contain a substantial fraction of intermediate-age stellar popu- 
lations and their prolon gued star formation hist ory is very similar 
to the picture drawn bv lAUard etal] j2006h and lSarzi etaU i2007h 




([S III] + Pa,)/Pa, 



Figure 7. Free emission line ratios versus photoionization models. The grid 
is computed, while the data points with eiTor bars represent the observa- 
tions. The metallicities and log(f7) ai'e identified. 

in the case of central star-forming rings, based on optical data. This 
similarity might support the view that central star formation often 
occurs in circumnuclear ri ngs, as was shown by the STIS Surve y 
of Nearby Nuclei (SUNNS ISarzi et al.ll2005LI Shields et al.ll2007l) . 
Analyses based on optical and NIR may be consistent in case of low 
dust obscuration or low contamination from old populations from 
the bulge region, which affect the optical colours. Our NIR based 
analysis is more robust in both respects, because of the whole sam- 
pling of the region between ~0.8 ^m and ~2.4 /im, which includes 
unique features of intermediate-age such as CN bands. 

We found that the metallicity of the stars which dominate the 
light at 1.2230/im, in the inner hundred parsecs of the four galax- 
ies is most likely solar. The reddening value determined with the 
whole spectra is lower than the one determined with the emission 
lines only. We found an excellent agreement of the free emission 
line spectrum with photoionization models, using as input SED 
the one derived via the spectral fitting in the NIR. Due to fea- 
tures like the l.l^mCN band, the stellar population analysis in 
the NIR seems to allow for a better fine tuning in the age dis- 
tribution of the galaxy stellar populations. The near-IR spectro- 
scopic approach is being now pursued a lso in other type of galaxies 
dSilva, Kuntschner & Lvubenovj|2008l) . 
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